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We have estimated the masses of the central supermassive black holes of 2442 radio galaxies 
from a catalog compiled using data from the NED, SDSS, and CATS databases. Mass estimates 
based on optical photometry and radio data are compared. Relationships between the mass of 
the central black hole M^^ and the redshift Zp are constructed for both wavelength ranges. The 
distribution of the galaxies in these diagrams and systematic e ects influencing estimation of the 
black-hole parameters are discussed. Upper-envelope cubic regression fits are obtained using the 
maximum estimates of the black-hole masses. The optical and radio upper envelopes show similar 
behavior, and have very similar peaks in position, Zp ~ 1.9 and amplitude, logMp''=9.4. This is 
consistent with a model in which the growth of the supermassive black holes is self-regulating, 
with this redshift corresponding to the epoch when the accretion-fiow phase begins to end and 
the nuclear activity falls off. 

PACS: 98.54.-h, 98.54.Gr, 98.62.Ve, 98.70.Dk, 98.80.Es 



1. INTRODUCTION 

One of the central questions in studies of distant ra- 
dio galaxies (redshifts z > 0.3) is the origin and 
growth rate of supermassive black holes (SMBHs) 
in their nuclei. The discovery of radio galaxies and 
quasars at high redshifts has raised the need to ex- 
plain the rapid formation of SMBHs in early epochs 
of the evolution of the Universe. For example, ra- 
dio galaxies have been detected at z = 5.19 (TN 
J0924-2201) T, when the age of the Universe was 
only t 1.1 X 10^ yrs (in a cosmology with il,n — 
0.27, nA = 0.73, H_0 = 70 km s-lMpc-1), and z 
— 4.515 [2j, when t ~ 1.3 x 10^ yrs. Quasars have 
been found in the Sloan Digital Sky Survey (SDSS) 
at redshifts of z = 6.42 (J114816.64-f 525150.3), 
z = 6.23 (J104845.05-h463718.3), and z = 6.05 
(J163033.90-f 401209.6) [3], when the age of the Uni- 
verse was t ^850-900 million years. Willot et al. [4] 
used the width of the emission lines of the quasar 
J114816.64-h525150.3 at z = 6.42 to estimate the 
mass of its SMBH to be Mbh = 3 x 10^ Mq. The 
rapid (over 800 million years) formation of SMBHs 
in a ACDM cosmology requires special models for 
their growth rates. As was noted by Loeb i_5^, if the 
formation of massive black holes is described by a 
heirarchical scheme similar to that for the formation 
of their parent galaxies, this process must be rapid 
and e cient. Models with primordial massive black 



holes that formed after the Big Bang but before the 
formation of galaxies are often invoked to explain the 
early existence of SMBHs. 

It is believed that SMBHs play a central role in the 
formation and evolution of massive galaxies, and are 
also a key component in the development of nuclear 
activity. However, it remains unclear how galaxies 
and their central SMBHs are related to the formation 
of observed structures [B]-[S]. An important result ob- 
tained from measurements of black-hole masses in 
nearby galaxies is the existence of an appreciable cor- 
relation between the mass of the central black hole 
and the mass of the bulge of the host galaxy, which 
gives rise to correlations between the black-hole mass 
and the luminosity of the bulge [5], [10] and between 
the black-hole mass and the stellar velocity disper- 
sion a (e.g. [ll]-[16]. This result is based on a small 
number of nearby galaxies for which direct black-hole 
mass measurements are available (^30), and there is 
appreciable scatter in the relation (^0.4 in the log- 
arithm of he black- hole mass). Nevertheless, these 
empirical relations provide a new tool for estimat- 
ing black-hole masses in various types of active galac- 
tic nuclei (AGNs), when the bulge luminosity and/or 
stellar- velocity distribution is known [T3-[IS]- The 
existence of these correlations can provide constraints 
on the evolution of AGNs, and can be used to test 
whether AGNs precisely follow these relations or not. 
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Franceschini et al. [2Dj established the existence of 
a correlation between the radio power of AGNs and 
the masses of their black holes. A number of studies 
have also demonstrated relationships between black- 
hole masses and radio powers for samples of nearby 
galaxies, but these results differ from one another 
[21] -[24]. Discussions of this relation were also ex- 
tended to more active galaxies, such as Seyfert galax- 
ies and quasars (see, for example, [^.[^-[25]). 

A number of studies have used data from the 
SDSS [29| for quasars to estimate the growth rates of 
their SMBHs. Kelly et al. [30, constructed the black- 
hole mass function for a sample of 9886 broadline 
quasars with redshifts 1 i z j 4.5. Their results support 
the "downsizing" effect for black holes in broad-line 
quasars, whose peak number density shifts toward 
high z as the black-hole mass is increasecQ. This peak 
occurs at z^2. Moreover, Kelly et al. [30] estimated 
the completeness of the SDSS sample as a function 
of black-hole mass and the ratio of the luminosity of 
a quasar to the Eddington luminosity, L/LEdd, and 
found that the sample of black-hole masses < 10^ 
and L/LEdd £ 0.5 is very incomplete at redshifts 
z I 1. According to their model, they estimate the 
lifetime of the broad-line quasar stage to be approx- 
imately 150 million years at z = 1 for a black- hole 
mass of ^ 10^ Mq. Shen et al. 31 constructed a sam- 
ple of 105783 quasars based on the SDSS, and used 
various properties (emission in Ha, H/3, Mgll, CIV; 
radio data; broad absorption lines) to estimate their 
black-hole masses. They also obtained a logarithmic 
dependence between the black-hole masses and the 
luminosities and widths of the lines used. 

Here, we verify two dependences of the black- hole 
masses on the radio and optical luminosities that are 
used in the literature in studies of radio galaxies. We 
used data from the SDSS, NED, and CATS databases 
to construct a sample of distant (z i 0.3) radio galax- 
ies, for which we derived and compared these depen- 
dences for the black-hole mass. We calculated the lu- 
minosities using optical data in the R filter and 5-GHz 
radio flux densities. We assumed a ACDM cosmology 
with HO = 71 km s-lMpc-1, 0^1/= 0.27, and fl^ = 
0.73. We used a catalog of distant radio galaxies with 
spectroscopic redshifts z i 0.3, which were selected 
based on data from the largest available databases 

[S2]-[3S]. 



The ''Downsizing'' scenario scenario, which can 
explain a mass dependence of the evolutionary history 
in which less massive elliptical galaxies have longer 
star-formation histories thain more massive elliptical 
galaxies . 



2. STUDIED SAMPLE AND MASS 
ESTIMATES 

2.1. Input Data 

A sample of radio galaxies with z i 0.3 
[32] -[35] was constructed using the NED 
(http : //nedwww . ipac . caltech. edu), CATS 
(http://cats.sao.ru) [35], [S3, and SDSS 
( |http : //www ■ sdss . org| ) [25] databases. This 
sample can be used to carry out various statistical 
and cosmological tests requiring a comparatively 
large number of objects of a single type [37]- [39]. 
The NED database was used to construct the 
initial list of sources, selecting objects with specified 
parameters, most importantly with redshifts z i 0.3 
and morphologies corresponding to radio galaxies. 
The initial catalog contained 3364 objects. This 
galaxy sample was contaminated by objects with 
incomplete information and objects with other 
properties. Therefore, we paid special attention to 
eliminating unwanted sources from the initial sam- 
ple, including (1) galaxies with redshifts determined 
using photometric methods and (2) galaxies with 
quasar-like properties according to data from the 
literature. The final catalog contains 2442 sources 
with spectroscopic redshifts, photometric data and 
radio flux densities, radio sizes, and radio spectral 
indices, calculated via cross identification with other 
radio catalogs in the CATS database at frequencies 
from 30 GHz to 325 MHz. 

2.2. Estimating the SMBH Masses. Optical 
Data 

To estimate the SMBH masss, we used the relation- 
ship between the black-hole masses and R absolute 
magnitudes {Mbh — A/r) from [16j . which was ob- 
tained via an analysis of the magnitudes and stellar 
velocities in AGNs: 

log{^) = -0.50(±0.02)AfR, - 2.27(±0.48) , 

Mq 

where Mr is the absolute R magnitude of the bulge, 
assuming that the bulge of an elliptical galaxy essen- 
tially represents the entire galaxy. 

Since R observations are not available for all the 
catalog radio galaxies, we estimated the R luminosi- 
ties of a number of objects based on data in other fil- 
ters using the database of galactic spectral-energy dis- 
tributions (SEDs) [40]. Photometric data were used 
to estimate the ages of the radio galaxies, assuming 
that their SEDs corresponded to the SEDs of ellipti- 
cal galaxies. We determined the R magnitude using 
the redshift and the age of the corresponding model 
SED track. We used the GISSEL model data (41j, 
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which contains synthetic spectra of elhptical galaxies 
for ages from 200 milhon years to 14 billion years, 
as the basic library of tracks. This library can be ac- 
cessed both in the HyperZ package [35] and at the 
site http : //sed . sao . ru| [40] . The distribution of es- 
timated and measured R magnitudes for the radio- 
galaxy sample is shown in Fig. 1. 

The luminosity was calculated using the formula 
[33] (see also [31) 

L = iirFdliz), 

z 

where F is the measured flux, dL{z) = {1 + z) J 



is the luminosity distance, and H{z) — 

7Jo[r!fl(i + + + z)3 - {no - i)(i + zf + rtA 

is the Hubble parameter. The resulting luminosity 
(more precisely, absolute magnitude) was used to 
estimate the masses of the central black holes of the 
galaxies. The distribution of estimated black-hole 
masses is shown in Fig. 1. 



2.3. Estimating the SMBH Masses. Radio 
Data 

To calculate the estimated black-hole masses based 
on radio data, we used the relationship between the 
radio luminosity and black-hole mass [P^ghz — Mth) 
from [20]: 

/MbftX 2.2-3.0 
^^^^^ " I Mo J 

where P5 qhz is the total power at 5 GHz. Our calcu- 
lations used the formula 

We calculated the radio power using the standard for- 
mula I43l 



P. 



4T:diS{ 



where s(^jj^:^^ is the source flux density at fre- 
quency ly taking into account the redshift and is 
the photometric distance. The value S'^-p^^^ was 
determined from the continuum radio spectrum of 
each source, constructed via fitting with a standard 
set of functions. We used the spg program [45] of 
the RATAN-600 continuum data-reduction system 
for this purpose. For radio galaxies with measure- 
ments at only one frequency, we took the spectral 
index to be equal to the median value for our sample 
[52] at 1.4 GHz {amed = -0.63). The distribution of 
the obtained radio power at 5 GHz is shown in Fig. 
2. 



3. DISCUSSION 

A comparison of the masses estimated using radio 
and optical data shows a clear difference in the shapes 
of the mass distributions (Figs. 1 and 2). Moreover, 
in contrast to the distribution of optical mass es- 
timates, the distribution of radio mass estimates is 
two-peaked. This could result from evolutionary man- 
ifestations in the radio flux at di erent redshifts at 
which the radio sources are ignited and then quiet 
down. Then, we may obtain di erent estimates for 
the mass of the "central engine" based on the ra- 
dio flux, corresponding to di erent stages of activity. 
The higher peak in the distribution in Fig. 2 corre- 
sponds to closer objects with lower masses, observed 

the SDSS. Note that, in spite of the two-peaked 
form of the mass distribution, the radio data provide 
a smaller scatter in the redshift distribution (Fig. 3). 
We should note that nonlinearity of "radio luminos- 
ity — optical luminosity" relation having the under- 
lying relation of radio luminosity and black hole mass 
was demonstrated in 1960-s by Pskovskii [4^ and cor- 
rected by Iskudarian and Pariiskii [47] . 

While the optical black-hole mass estimates can 
tentatively be divided into two subgroups-strongly 
clustered at z j 0.7 and weakly clustered at z i 0.7- 
this is not the case for the radio mass estimates. The 
range of black-hole masses derived using the radio 
data is two orders of magnitude smaller than that for 
the optically derived masses. The mean radio masses 
in individual redshift intervals {Az = 0.5) are higher 
than the corresponding optical masses. Although the 
radio estimates also display clustering at z j 0.7 for the 
flrst subsample, the objects of the second subsample 
are grouped along the line for the median regression 
fit. 

Figure 3 shows the least-squares regression fits 
and upper envelopes for both types of data. For the 
optical data (Fig. 3, left), the fit for the median val- 
ues on the log MRbh-z diagra in bins Az = 0.5 
yields a growing function satisfying the linear rela- 
tion logMf = 0.387(±0.134) • z + 7.665(±0.309), 
where is the black-hole mass estimated using 

the R optical data. The upper envelope is described 
by a cubic function: logM^'' = 9.058 + 0.783 • z - 
0.303 • + 0.030 • z^. Similarly, the linear fit for the 
radio data (Fig. 3, right) has the form logM^^^ — 
0.126(±0.027) •z-H8.793(±0.073), the upper envelope 
the form logM^^^^ = 9.084 + 0.347 • z - 0.120 • z'^ + 
0.011 • z^. In spite of the difference in the dispersion 
of the mass estimates, the posi-tions and amplitudes 
of the maxima of both upper envelopes are similar: 



Here and below, the black-hole masses are expressed 
in solar masses, M0 
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Figure 1: Left: distribution of R magnitudes for the sample of radio galaxies. Right: distribution of black-hole 
masses obtained from R data for the sample of radio galaxies. The SDSS subsample of radio galaxies is shown 
by the dotted line. 
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Figure 2: Left: distribution of the 5-GHz radio power for the radio-galaxy sample. Right: distribution of the 
black-hole masses obtained from the total radio power (luminosity) at 5 GHz. The SDSS subsample of radio 
galaxies is shown by the dotted line. 




the peak is at Zp = 1.78 and logM"^ = 9.67 for the 
optical data, and at Zp — 1.92 and log A/p'*=9.38 for 
the radio data. 

The location of the maxima in the redshift range 
1.5 i z i 2 (with similar values for both the optical and 
radio estimates) could be associated with some sort 
of systematic e ect reflecting real physical pro-cesses. 
This range of z corresponds to an epoch of massive 
"ignition" of radio sources as a result of mergers be- 
tween galaxies in galaxy clusters. In this case, due to 
selection effects, we will most likely detect objects 
with the maximum luminosity, which means those 
with the maximum estimated black-hole masses, in 
this range. 

Note that this result obtained for our list of radio 
galaxies is consistent with the results of [30] , where a 



sample of broad-line quasars from the SDSS was con- 
sidered. Kelly et al. [30] found that the peak num- 
ber density of SMBHs in broad-line quasars occurs 
at z 2. The location of the mass-distribution peak 
at z ~ 2 is consistent with a self-regulating growth 
model for the SMBHs, with this being the epoch when 
the end of the accretion-flow phase begins and the nu- 
clear activity starts to fall off. Note that our estimates 
and the estimate of [30] are also consistent with self- 
regulation of the rate of formation of galaxies in halos 
of dark matter via the interaction of the jets from the 
central object with the intergalactic gas [48] . 

A growth in the black-hole mass with increasing 
redshift is visible in both diagrams, demonstrated by 
the linear fits. The uncertainties in the fit parame- 
ters characterize the statistical scatter in the masses 
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Figure 3: Plots of black-hole mass versus redshift for the radio galaxies, based on R optical data (left) and 
5-GHz radio data (right). Regression fits were constructed for the mean and maximum black-hole masses in 
intervals Az = 0.5. 



(and luminosities) derived using the optical and radio 
data. The large scatter for the optical estimates, in 
particular the many deviations from the maximum 
values, indicates small ratios of the optical to the 
Eddington luminosities, and may suggest that some 
Seyfert 2 galaxies have been included in our sample. 
The linear fit based on the median masses may re- 
flect selection effects in the detection of distant ob- 
jects, and the region of intersection with the upper 
envelope (zc = 4.14,log(Mc) = 9.27 for the optical 
and Zc = 3.57,log(Mc) = 9.25 for the radio) shows 
an upper limit, beyond which this effect will domi- 
nate. Note, however, that only a few radio galaxies 
have been detected in this region |2]. 



The M„^p^ 



rad 



relation indicates the region in 
which the two methods used to estimate the masses 
of the SMBHs are statistically applicable, and shows 
the amount of scatter obtained using these methods 
(Fig. 4). Three regions of clustering can be distin- 
guished on this plot, for two of which we obtained 
fits in the form log M^j^ = a ■ log M^^^^ + b The pa- 
rameters for the first group of clustered points are 
a = -0.036 ± 0.026, b = 8.122 ± 0.219, and those 
for the second group of points a = 1.211 ± 0.029, 
b = —2.329 ± 0.261. The data for radio galaxies from 
different redshift ranges are shown by different sym- 
bols: 0.3 < z < 0.7 (points), 0.7 < z < 1.5, (crosses), 
and 1.5 < z (triangles). The high concentration of 
points in the left part of the diagram is due to the 
large number of comparatively nearby radio galaxies 
in the SDSS, for which relatively similar black- hole 
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Figure 4: Plot of M^^t versus M^^^ for R-band and 
5- GHz data. Regression fits were obtained for two 
regions where the points are concentrated. Data for 
radio galaxies from di erent redshift ranges are de- 
noted with different symbols: 0.3 < z < Q.7 (points), 
0.7 < z < 1.5 (crosses), and 1.5 < z (triangles) . 



masses were obtained. 

Figure 4, and also Figure 1, show that the mean 
estimated masses for the "central engines" for the 
SDSS subsample of galaxies are ^ 10* Mq based on 
the optical data and from 10* to 10^ Mq based on the 
radio data. In general, nearly all the radio galaxies in 



the sample have M^^^ 



— ^rad- 1^0^ ^ second group. 



many of which are distant radio galaxies with z > 1.5 
we observe a correlation between the two mass esti- 
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mates for the central objects. The objects in the fower 
right part of Fig. 4 have fow optical and high radio 
mass estimates. These low estimates from the region 
of the lower envelope on the left Mbh(z) diagram in 
Fig. 3 and the corresponding estimates in the right 
diagram may be associated with non-elliptical sub- 
classes of radio galaxies. 

We can distinguish several selection e ects that 
are inevitable when estimating the central black-hole 
masses, even when care is taken in selecting galax- 
ies according to their optical and radio properties. 
Their manifestation can explain differences between 
results obtained using different methods to determine 
the central black- hole masses. These selection effects 
include the following: 

1) more distant objects tend to be more powerful; 
i.e., a sample of distant objects contains, on average, 
galaxies with higher black-hole masses; 

2) although the sample has been cleaned of ob- 
jects with the properties of quasars (such as broad 
Ha emission), some such objects with relatively weak 
lines may remain in the subsample of comparatively 
nearby (z j 0.7) sources; 

3) the SDSS objects in our catalog should not in- 
clude Seyfert 2 galaxies with narrow lines; in spite of 
their radio emission, these galaxies have low ratios of 
their bolometric to Eddington luminosities (see, e.g., 
[49]), and their luminosities should thus place them 
in the lower part of the diagram (Fig. 2); 

4) the presence of evolutionary effects that are 
strongly manifest in the radio, such as the fact that 
some radio galaxies may be observed in early phases 
of their activity, and that galaxies with central black 
holes of the same mass may have radio fluxes differing 
by an order of magnitude; 

5) some comparatively nearby objects whose peak 
activity in the radio has passed have black-hole mass 
estimates derived from velocity dispersions calibrated 
based on the radio fluxes; in this case, the use of de- 
pendences for more distant galaxies can lead to over- 
estimation of the black-hole masses. 

Nevertheless, the similarity of the behaviors of the 
mass upper envelopes (Fig. 3) suggests that these 
estimates are trustworthy for objects with the high- 
est luminosities (close to the Eddington values). The 
comparatively small scatter of the radio masses for 
various z makes it possible to use this diagram to es- 
timate masses with a mean accuracy of less than an 
order of magnitude, or even half an order of magni- 
tude at z I 1.5. This is confirmed by the plot in Fig. 
4, where a correlation between the radio and optical 
masses is observed for most radio galaxies with z i 
z > 1.5. 

We used the other parameters in our catalog to 
construct plots of the black-hole mass versus radio 



spectral index (Fig. 5) and versus the apparent size 
of the radio source (Fig. 6) for a single type of galax- 
ies with active nuclei (i.e., radio galaxies). The former 
plot shows a correlation between the radio mass esti- 
mates and the radio spectral index 

a = 2.73(±1.56) - 0.18(±0.08) • log{Mbh) , 

There is a tendency for more powerful radio galax- 
ies, which presumably have higher nuclear activiity, 
and possibly higher masses for their central engines, 
to have steeper spectra (i.e., smaller spectral indices). 
Of course, it cannot be ruled out that some selection e 
ect is hidden via the redshift dependence of the spec- 
tral index a{z) found in [35j. Distant radio galaxies 
are often selected in surveys based on their spectral 
indices, and distant sources tend to be more powerful; 
therefore, this e ect should make some contribution to 
the detected dependence. However, the physical prop- 
erties of the objects themselves at various z probably 
dominate in the a{z) relation. 

2.5 I — ■ ■ ■ ■ — ■ ■ ■ 1 

1.5 - ' . 




-^.3 ' ' ' — ' — ^ ' ' ' ' ' ' — ' — "5 ' ' ' ' 

10 10 

M„"' 

Figure 5: Plot of black-hole mass versus radio spec- 
tral index. The masses were estimated using the radio 
data. Notation is the same as in Fig. 4- 

Figure 6 shows the distribution of the angular 
sizes of radio galaxies according to the NVSS data [SO] 
as a function of their radio black-hole mass estimates. 
Most of the objects distributed along the horizontal 
axis have small angular sizes (i25"), determined by 
the resolution of the NVSS (~45"), and most of these 
are SDSS objects with redshifts z j 0.7. Note that the 
central objects of the most extended objects {l 1') 
cover the entire range of masses, and are encountered 
in all three redshift intervals we have considered. 

4. CONCLUSION 

We have carried out a comparative analysis of esti- 
mates of the central black-hole masses of 2442 radio 
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Figure 6: Masses estimated using the radio data. The 
apparent sizes of the galaxies were taken from the 
NVSS catalog fj^. Notation is the same as in Fig. 

4- 

galaxies with z i 0.3 [32]-|33]: derived from relations 
between the black-hole mass and the R luminosity 
[T7] and between the black-hole mass and the radio 
power [22] ■ Appreciable differences between these two 
estimates are observed for many of the radio galaxies. 
We have discussed various systematic effects associ- 
ated with observational selection effects and the evo- 
lutionary properties of the radio galaxies that could 
lead to these differences. However, a diagram of AI^ 
versus M^^^ (Fig. 4) reveals a region where these two 
mass estimates are correlated. This zone is formed 
primarily by the distant radio galaxies in our sam- 
ple. Moreover, the upper envelopes constructed us- 
ing the maxima of the two mass estimates (Fig. 3) 
show similar behavior and have very similar posi- 
tions {zp ~ 1.9) and amplitudes (log M^''=9.4). This 
is consistent with self-regulation of the growth of the 
central SMBH, when the end of of the accretion flow 
begins and the activity of the galactic nucleus falls off. 
The M^^^(z) diagram displays comparatively narrow 
scatter, and should be preferred for use in estimating 
galactic black-hole masses. 
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